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ABSTRACT 

We study the internal dynamical structure of the early-type lens galaxy SDSS J0728+3835 
at z = 0.206. The analysis is based on two-dimensional kinematic maps extending out to 
1 .7 effective radii obtained from Keck spectroscopy, on lensing geometry and on stellar mass 
estimates obtained from multiband Hubble Space Telescope imaging. The data are modelled 
under the assumptions of axial symmetry supported by a two-integral distribution function 
(DF), by applying the combined gravitational lensing and stellar dynamics code cauldron, 
and yielding high-quality constraints for an early-type galaxy at cosmological redshifts. Mod- 
elling the total density profile as a power-law of the form ptot , we find that it is nearly 
isothermal (logarithmic slope y' = 2.08^qq2)' ™d quite flattened (axial ratio q = 0.60^q q3). 
The galaxy is mildly anisotropic (6 = 0.08 + 0.02) and shows a fair amount of rotational sup- 
port, in particular towards the outer regions. We determine a dark matter fraction lower limit 
of 28 per cent within the effective radius. The stellar contribution to the total mass distribu- 
tion is close to maximal for a Chabrier initial mass function (IMF), whereas for a Salpeter 
IMF the stellar mass exceeds the total mass within the galaxy inner regions. We find that the 
combination of a NFW dark matter halo with the maximally rescaled luminous profile pro- 
vides a remarkably good fit to the total mass distribution over a broad radial range. Our results 
confirm and expand the findings of the SLACS survey for early-type galaxies of comparable 
velocity dispersion (itsdss = 214 + 11 km s"'). The internal structure of SDSS J0728 is con- 
sistent with that of local early-type galaxies of comparable velocity dispersion as measured 
by the SAURON project, suggesting lack of evolution in the past two billion years. 



Key words: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: kine- 
matics and dynamics — galaxies: structure. 



1 INTRODUCTION 

Unveiling the mass distribution, dynamical structure and dark mat- 
ter content of early-type galaxies is of great interest both as a sub- 
ject in its own right, considering their importance in the local Uni- 
verse, where a large fraction of the total s tellar mass is contained 
within E/SO systems jpukugita. Hogan."&^ Peebles 1998), and in or- 
der to provide stringent tests for the galaxy formation and evolution 
models. 

It is not surprising, therefore, that nearby early-type galax- 
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ies have been the object of intense study during the last decades, 
by taking advantage of the diverse avail able observational tracers. 



These include stellar kinematics (see e.g.lSaglia. Bertin. & Stiavelli 



^agl 

1992', 'Franx et al.' 'l994', 'Rix e t all 1 19971 'Loewenstein & White 
1999, Gerhard etal 2001 Borriello et a l j ,2003. Cappellari et ~ 



20071 iThomas et al.ll2007bl . 1weiimans et al.H2009l), globular clus_ 
ters and planetary nebul a e kinematics (e. g . Icote et al] l2003l : 
iRomanowskv et aU l2003l : Ide Lorenzi et al.l l2008h. the occa- 
siona l HI disk or ri n g (e.g . Franx. van Gorkom. & de Zeeuwl 
1994 IWeiimans et all \ 200^ 



^ and hot X-ra y emission (e.g. 

Matsushita et alj Il99 8: Fukaz awa et alj l2006l : iHumphrev et all 
20061: iHumphrev & Buota |2010|) . The general picture emerging 
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from many of these studies is that the total mass density profile 
of ellipticals can be well described by a power-law form close to 
Plot l/r", often referred to as the isothermal profile. Moreover, 
while the inner regions of early-type galaxies are clearly dominated 
by the stellar component, the dark matter component is usually 
found to play already a non negligible role, with fractions of ap- 
proximately 10 to 40 per cent of the total mass within an effective 
radi us. Studies based on stellar population and dynamical models 
(e.g. IPadmanabhan et al]|2004) indicate that the dark matter frac- 
tion increases with the mass of the galaxy, a tr end that is more con - 
spicuous in the case of slow-rotator ellipticals l lTortora et al.l2009l) . 

The analysis of early-type galaxies beyond the local Universe, 
i.e. beyond redshift z ~ 0.1, holds great promise in view of un- 
derstanding the structural evolution of these objects, but it also 
presents several difficulties which hinder the application of tradi- 
tional techniques. Stellar dynamics studies, in particular, are lim- 
ited by the degeneracy between the galaxy mass profile and the 
anisotropy of the stellar velocity dispersion tensor. Taking into 
acc ount higher order velocity moments can provide a solution 



(seelGerhar3ll993l : Ivan der Marel & Franxl 19931 : iLokas & MamonI 
l2003h . but unfortunately carrying out such measurement for distant 
systems is not viable with the current instruments on 8-10 meter 
class telescopes. 

However, galaxies at z > 0.1 have a far greater chance of act - 
ing as strong gravitational lenses ( [Turner. Ostriker. & GottllT984l) . 
thus providing a very helpful additional diagnostic tool. This is 
particularly valuable since it allows an accurate and robust de- 
termination of the total mass enclosed, in projection, with in the 
region delimited by the Einstein radius l lKochanekl [T99ll) . Un- 
fortunately, the diagnostic power of strong lensing to constrain 
internal mass distribution of the deflector is l imited, chiefly by 
the mass-sheet and mass-slope degeneracies jpalco et al] 1 19851 : 
IWucknit j |2002|) . although the latter can be part ially overcome 
by st udying spatially extended lensed sources (e.g. I Warren & Dv^ 
12003 : Suyu et al. 2010). A very effective way to overcome 
these difficulties and robustly recover various structural proper- 
ties of the galaxy is to combine the gravitational lensing anal- 
ysis with the complementary constraints pro v ided by stellar dy- 
namics (see Koopmans & Treu 2002, 2003, Tre u & Koopmana 
[2002a,200 2b, 2,004, Barnabe & Koop m gns 2007, hereafter BK07 
and, e.g.. iRusin & Kochanekl I2005L Jjiang & Kochanekl l2007l 



Ivan de Ven et al.ll2008l , iTrott et alJl20ialGrillo et alj|2010l for fur- 
ther applications of this approach). 

Up to very recently, the availability of only a handful of lens 
galaxies suitable for the joint analysis represented a major limi- 
tation. This con cern has been dispelled by the Sloan Lens ACS 
Survey, SLAC S (iBolton et al.'2006V2008a'.'2008bVKoop mans et al .1 
2006 , 2009 . Treu et alj 12006. 2009, Gavazzi et al. 20071 |2008| . 
Auger et al.l2009l) , which has led to the discovery of a large and ho- 



mogeneous sample of almost a hundred strong gravitational lenses, 
mostly early-type galaxies, in the redshift range of z ~ 0.05 - 0.5. 
For a subset of about 30 SLACS systems, the data set is com- 
plemented by two-dimensional kinematic maps of the lens ob- 
tained from spectroscopic observations carried out either with the 
Very Large Telescope (VLT) instrument VIMOS or with the Low 
Resolution Imager and Spectrograph (LRIS, see lOke et alj[T993) 
mounted on the Keck-I telescope. This has provided further mo- 
tivation to expand the combined analysis technique into a more 
general and self-consistent method which makes full use of the 
available data sets (i.e. surface brightness distribution of both the 
lensed source and the lens galaxy, and two-dimensional kinematic 
maps of the latter), and is coherently embedded in the framework 



of Bayesian statistics (BK07). The current implementation of the 
method, the cauldron code — based on the assumptions of axial 
symmetry and two-integral stellar DP for the lens galaxy — has 
been used to conduct an in-depth study of a SLACS subsample 
of six systems representativ e of the survey in ter ms of redshifts 
and ve locity dispersions (see lCzoske et alj[2008l and lBamabe et al.l 
l2009al hereafter B09). As shown in those works, the more sophisti- 
cated approach makes it possible to extract much more information 
out of the data set, allowing to recover, in addition to the slope of 
the total density profile, several other important properties of the 
lens galaxies, including the flattening of the density distribution, 
lower limits for the dark matter fraction at different radii and in- 
sights on the dynamical structure (angular momentum, anisotropy, 
contribution of rotation and random motions). These quantities are 
all of relevance to the formation history of these galaxies. 

In this paper we carry out a detailed combined lensing and 
dynamics analysis of the SLACS system SDSS J0728-I-3835, em- 
ploying the CAULDRON algorithm. The lens is an early-type galaxy 
at z = 0.206, with an aperture averaged velocity dispersion a = 
214 ±11 km s"' measured from SDSS spectroscopy and a half- 
light radius = 1.78" in the / band. The background source is 
located at z = 0.688 and the Einstein r adius is i?Eins i = 1 .25". 
With respect to the systems considered in lCzoske et al. ( |2008|) and 
B09, all followed-up with VLT VIMOS integral-field unit, the main 
difference in the observables lies in the kinematic data set: the 
SDSS J0728 velocity moments maps are obtained from LRIS Keck 
long-slit spectroscopic observations, using three slits parallel to 
the major axis, and offset along the minor axis, in order to mimic 
integral-field capabilities. Remarkably, this kinematic data set ex- 
tends significantly farther than those of the previously examined 
systems, reaching up to 1.7 R^, thus providing us important con- 
straints beyond the inner regions of the galaxy. Moreover, for the 
first time, we us e the stellar masse s determined from stellar pop- 
ulation analysis jAuger et al.ll2009l) to set the normalization of the 
luminous mass distribution of the lens galaxy, enabling us to disen- 
tangle the luminous and dark matter contributions and to compare 
different choices of the initial mass function (IMF). 

The paper is organized as follows: after introducing the data 
set in Section|2l we present and discuss the results of the combined 
analysis in Section|3]and draw conclusions in Section|4] 

Throughout this paper we adopt a concordance ACDM model 
described by Dm = 0.3, Q.^ = 0.7 and tfo = lOO/zkms"' Mpc"' 
with h = 0.7, unless stated otherwise. 



2 OBSERVATIONS 

2.1 High-resolution imaging data 

The lensing analysis requires deep high-resolution imaging data 
and this is provided by the SLACS survey's Hubble Space Tele- 
scope (HST) imaging. In particular, SLACS has obtained one HST 
orbit (a 2200s) of d ata in the P814W fi lter. The data are pro- 
cessed as described in iBolton et all ( l2008al) : to briefly summarize, 
the four individual exposures are background subtracted, cosmic 
ray-cleaned, registered, resampled to an output grid with square 
pixels that are 0'.'05 on a side, and stacked with an additional cos- 
mic ray-rejection step. Synthetic point-spread function (PSP) im- 
ages created with Tiny Tim are likewise resampled and combined 
to create a composite model PSP for the output image. 

The li ght distribution of th e galaxy is then fit with a B-spline 
model (e.g. . IBolton et al .l2008al) and this model, convolved with the 
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Figure 1. HST F814W image of SDSS J0728 with tlie slit locations over- 
plotted as solid red lines. The most narrow slit ((X'V) is centered on the light 
distribution and follows the semi-major axis. The other slits are l"and I'.'S 
wide and are offset so as to be adjacent to the previous slit. The dotted red 
lines delineate the 7 apertures with width (XHS that were extracted from 
each longslit observation. 



with a C.'T slit positioned on the center of the galaxy; four 1200s 
exposures with a l"sht offset 0'.'85 from the center of the galaxy; 
and 4 1800s exposures with a I'.'S slit offset 2'.'1 from the center 
of the galaxy. The data were re duced using custom Python scripts 
(for details see ISuvu et alj2o"Tol) and one dimensional spectra were 
extracted from apertures of width 2 pixels (« 0'.'43) at seven points 
along each slit as indicated in Figure [T] yielding 21 spectra in total 
(e.g., Figure|2ll. 

The stellar velocity disper sion and ve l ocity offset was com- 
puted for each spectrum as in ISuvu et al.l ( 1201 Ol) using the rest- 
frame wavelength range 4200 A to 5450 A and errors were deter- 
mined from Markov chain Monte Carlo simulations. The resulting 
velocity and velocity dispersion profiles for each slit position are 
shown in Figure [3] We compare our data with the velocity disper- 
sion derived from the SDSS spectrum by combining the 14 spectra 
from the inner two slits and determining a composite velocity dis- 
persion. This is found to be 236 km s"', rather larger than the value 
of 214 km s"' found for the SDSS spectrum. We have re-analysed 
the SDSS spectrum and find a = 212 km s"', but if we mask the 
O II lines (as is done in the analysis of the Keck spectra) we find 
cr = 225 km s"' . This is still somewhat lower than, although consis- 
tent within the errors with, the central velocity dispersion derived 
from the composite Keck spectrum. However, we might expect the 
SDSS value to be lower due to poor seeing pushing more flux from 
large radii (where the velocity dispersion is lower) into the SDSS 
fibre aperture. 



seeing present during the spectroscopic observations, is used as the 
input surface brightness distribution for the dynamics modelling. 
The residual image (the B-spline model subtracted from the data 
image) contains the flux from the lensed background source and 
is used as lensing constraint on the mass model of the lens galaxy 
(see BK07 for details). The galaxy-subtracted image is shown in 
the bottom right-hand panel of Fig.|4] 

2.2 Pseudo-IFU Spectroscopy 

Our previous joint lensing and dynamics studies used VLT-VIMOS 
integral field unit (IFU) spectroscopy to generate resolved stel- 
lar kinematic maps. In this study we employ 'pseudo-IFU' spec- 
troscopy, which consists of dithering a traditional longslit in the 
spectral direction over the target galaxy to provide spatially- 
resolved kinematic data perpendicular to the slit direction in ad- 
dition to the spatial information that is obtained along each slit (see 
FigurelTJ. This technique yields significantly fewer spaxels per ex- 
posure than proper IFU spectroscopy but the throughput from the 
longslit is a factor of a few larger than the IFU throughput and we 
are therefore able to efficiently produce kinematic maps, despite the 
less efficient way of sampling the kinematic field. Additionally, we 
are more sensitive to the low-surface brightness outer regions of the 
galaxy and can therefore extend our kinematic maps to larger radii 
than is possible with the IFU data. 

The spectroscopic data for SDSS J0728 were obtained with 
LRIS on Keck I during the nights of 22 and 23 December 2006. 
The observing conditions were clear with w O'.'B seeing, and we 
used the 460 dichroic to split the beam to the blue and red sides 
of the spectrograph. Here we only use data from the red side, 
which employed the 900/5500 grating with a dispersion scale of 
0.85 Apixel"'. Three slit positions aligned along the semi-major 
axis of the galaxy were observed, including: two 900s exposures 



3 ANALYSIS AND RESULTS 

3.1 Joint gravitational lensing and stellar dynamics analysis 

Here we briefly recall the general features of the cauldron code, 
the tool employed to carry out the combined self-consistent gravita- 
tional lensing and stellar dynamics analysis for the galaxy in exam, 
SDSS J0728. The reader is referred to BK07 for a detailed descrip- 
tion of the algorithm. 

We characterize the lens galaxy by means of its total density 
distribution p(//), where // is a set of parameters describing the den- 
sity profile. Via the Poisson equation, we calculate the total gravi- 
tational potential O and we use it simultaneously for both the grav- 
itational lensing and the stellar dynamics modelling of the data set, 
which typically includes the surface brightness and velocity mo- 
ments maps of the lens galaxy and the surface brightness map of 
the lensed image. Both these modelling problems can be formal- 
ized as a set of regularized linear equations, for which standard 
solving techniques are readily available. Thus — given a combined 
data set — for each choice of the parameters i] we can calculate 
the surface brightness distribution of the unlensed source, and the 
weights of the elementary stellar dynamics building blocks. In or- 
der to determine the "best" (in an Occam's razor sense) density 
model given the data, this linear optimization scheme has been em- 
bedded within the framework of Bayesian statistics. This allows to 
objectively quantify (and therefore rank) the plausibilit y of each 
model by m eans of the evidence merit function (see e.g. iMacKavl 
Il999ll2003h . In this way, by maximizing the evidence, one recovers 
the set of non-linear parameters // corresponding to the best density 
model, i.e. the model which maximizes the joint posterior proba- 
bility density function (PDF), hence called maximum a posteriori 
(MAP) model. 

The method as described is extremely flexible and can in 
principle support any density profile, by adopting for example a 
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Figure 2. Extracted one-dimensional spectra for the central aperture of the (Y.'l slit (left), l"slit (center), and I'.'S slit (right). The red line is the best-fit model, 
the green fine is the polynomial continuum model, and the lower panel shows the residuals. Greyed-out regions indicate parts of the spectra that were not 
included in the fit (two strong night-sky features, the Mgb line of SDSS J0728, and the O II features from the background source). 
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Figure 3. Line-of-sight velocity (left) and velocity dispersion (right) profiles for the lens SDSS J0728. The black circles are for the central slit, the blue squares 
are for the slit offset ff.'85 from the center, and the red diamonds are for the widest slit offset 2'.' 1 from the galaxy center. 



completely general pixelized density distribution. However, its cur- 
rent practical implementation, the cauldron algorithm, is more re- 
stricted in order to make it computationally efficient and assumes 
axial symmetry (i.e. a density distribution of the form p(R,z)) 
and a two-integral stellar DF / = f{E, L-), where E and L- de- 
note the two classical integrals of motion, i.e., respectively, en- 
ergy and angular momentum along the rotation axis. As shown 
in BK07, under these assumptions it is possible to take advan- 
tage of a fast Monte Carlo numerical imple mentation of the two- 
inte gral Schwarzschild method described by ICretton et al.l j 19991) 



and lVerolme & de Zeeuwl j2002h . which allows a dynamical model 
to be built in a matter of seconds and, therefore, makes it possible 
to explore large regions of the parameter space. The distinguishing 
feature of this method is that the building blocks employed for the 
construction of the dynamical model are constituted by two-integral 
components (TICs) ra ther than stellar orb its as i n the classical 
Schwa rzs child method i Schwarzschil j 1 979l ; see e.g. lThomas et al.l 
l2007bl and lvan den Bosch et al. l2008l for modem implementations). 
A TIC can be visualized as an elementary toroidal system, com- 
pletely specified by a particular choice of E and L-. TICs have sim- 
ple 1/R radial density distributions and analytic unprojected veloc- 
ity moments, and by superposing them one can build f(E, L,) m od- 
els for arbitrary spheroidal potentials (cf. ICretton et alj[l999l) : all 
these characteristics contribute to make TICs particularly valuable 
and inexpensive building blocks when compared to orbits. 



Notwithstanding these restrictions, iBamabe et alj ( l2009lj) 
have shown that cauldron works robustly even when applied to 
simulated systems which depart significantly from the method's as- 
sumptions (including the assumption of axial symmetry), reliably 
recovering several important global properties of such systems, in 
particular the slope of the total mass density profile, which is de- 
termined within less than 10 per cent of the correct value. When 
the system displays rotation in the kinematical maps, as is the case 
for SDSS J0728, the main dynamical quantities (such as the global 
anisotropy parameter <5, the angular momentum and the ordered to 
random motions ratio, see Sect. 13. 5t are recovered with an accuracy 
of 10 to 25 per cent. 



3.2 The galaxy model 



Gerhard et alj|200ll). s t rong and weak gray 



Koopmans et al.ll2009l: lOavazzi et al]|2007h 



Stellar dynamics (e.g. 

itational lensing (e.g. =^ 

and X-ray studies (see e.g. Humphrev & Buot3 20 id . and refer- 
ences therein) all concur in indicating that the total mass profile of 
elliptical galaxies is remarkably well described by a single power- 
law model over a large radial range. 

As our model for the total mass density profile of the analyzed 
lens galaxy we adopt, therefore, an axially symmetric power-law 
distribution 
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p(m) : 



my' 



(1) 



with po being a density scale, < y' < 3 the logarithmic slope of 
the density profile, and m the elliptical radius, i.e. 



Table 1. Recovered model parameters for lens galaxy 
SDSSJ0728. 



R 



2 +3 



(2) 



where Cg and ag are length-scales and the axial ratio q = Cg/ao. 
The corresponding (inner) gravitational potential (^(R,z) associ- 
ated with the distribution p(m) can be obtained straightforwardly 
by means of the classical IChandrasekhaJ i ll 9691) formula, which 
entails the computation of a single integral. 

In case the assumption of a power-law density profile is 
an oversimplified description of the true mass distribution of the 
galaxy, this is expected to have visibly disrupting effects on the 
reconstructed ob servables, in particula r for the lensing ones. As 
illustrated in the iSamabe et al. I l l2009bh tests, these can include a 
strongly irregular reconstructed source, with sharp transitions in 
intensity between close pixels (despite the usual adoption of a reg- 
ularization term), and the presence of recognizable arc-like fea- 
tures in the lens image residuals. We emphasize, however, that such 
clear discrepanci es have never surfac ed in previous analyses of the 
SLACS systems jCzoske et ai]|2008l B09), which supports instead 
the effectiveness of the simple power-law model. 

With this choice, the galaxy model is therefore characterized 
by three non-linear physical parameters, i.e. the slope y', the axial 
ratio q and the dimensionless lens strength Uq, which is directly re- 
lated to the normalization of the three dimensional potential (see 
Appendix B of BK07). To these, we must add the four "geometri- 
cal" parameters defining the configuration of the system in the sky, 
i.e. the position angle j?pa, the inclination / and the coordinates of 
the lens galaxy centre. The latter, as well as the angle i?pa, are typi- 
cally strongly constrained by the lens image brightness distribution 
and can be accurately determined by means of fast preliminary ex- 
plorations and, therefore, kept fixed afterwards in order to decrease 
the number of free non-linear parameters during the more compu- 
tationally expensive optimization and error analysis runs. If neces- 
sary, external shear can be also accounted for, by introducing shear 
strength and shear angle as additional parameters. 

Finally, we have three so-called hyperparameters which con- 
trol the level of the regularization in the reconstructed quantities: 
one for the surface brightness distribution of the unlensed source, 
and two for the TIC weights map. Their values are set by optimiz- 
ing the Bayesian posterior probability. 

3.3 Uncertainties 

In order to correctly assess the model uncertainties within the 
framework of Bayesian statistics, one needs to evaluate the pos- 
terior probability distribution of the parameters, i.e., by denoting 
the data set as d and the considered model or hypothesis as H, 



nr]) = Mri\d,'H)oc£(n)xp(rj), 



(3) 



where £(i]) = Pr(rf|j7,'K) is the likelihood and p(i]) = Pr{i]\'H) 
is the prior. The individual parameter confidence intervals can be 
obtained by marginalizing the posterior over each of them. 

Sampling a multidimensional distribution such as P is in gen- 
eral a challenging and computationally expensive task. An ef- 
fective technique to tackle this probl em is th e neste d sampling 
Monte Carlo method introduced by ISkillind ( |2004|) which, in 
calculating the evidence, produces posterior inferences as valu- 
able by-products. For our error analysis, we make use of the 
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2.082 


2.055 


2.077 


2.119 




0.325 


0.317 


0.323 


0.331 
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0.602 


0.574 


0.603 


0.681 


i 


68.1 


67.8 


68.6 


73.6 



Note. The listed parameters are: the logaiithmic slope y': 
the dimensionless lens strength ag; the axial ratio q and 
the inclination (' (in degrees) with respect to the line-of- 
sight. The second column presents the MAP parameters, i.e. 
the parameters that maximize the joint posterior distribu- 
tion. Columns 3 to 5 encapsulate a description of the one- 
dimensional marginalized posterior distributions, plotted in 
Fig. |6] The parameter ;/,^ax corresponding to the maximum 
of that distribution is listed in column 4, while and 
ai'e, respectively, the lower and upper limits of the 99 per 
cent confidence interval. 



MultiNest algorithm developed by iFeroz & HobsorJ l l2008t) and 
iFeroz. Hobson. & Bridges! dioO^ , which provides an efficient and 
robust implementation of the nested sampling method, and has been 
shown to yield reliable posterior inferences even in presence of 
multi-modal and degenerate multivariate distributions. 

The model parameters that we consider are the ones in- 
troduced in the previous Section, i.e. the inclination /, the lens 
strength ag, the slope y' and the axial ratio q (the additional param- 
eters which can be estimated by means of preliminary runs, such as 
the lens center and the position angle, however, are kept fixed here 
in order to reduce the computational load), to which we must add 
the three hyper-parameters. We formalize our ignorance by adopt- 
ing flat priors (or flat in logarithm for the hyper-parameters), con- 
structed around the MAP value of each parameter, and taken wide 
enough to include the bulk of posterior probability distribution. Fi- 
nally, we conduct a thorough exploration of this 7-dimensional sur- 
face by launching MultiNest with 400 live pointfl from which 
the individual marginalized posterior probability distributions are 
obtained. In the following, we quote the 99 per cent confidence in- 
terval calculated from these distributions as our error. 



3.4 Results for the best model reconstruction 

The CAULDRON code has been applied to the combined data set de- 
scribed in Section |2] In order to avoid possible spurious effects in 
the lensing reconstruction, we have masked out the central regions 



The live points (sometimes also called active points) are the initial sam- 
ples drawn from the full prior distribution p(i]), from which the nested 
sampling exploration of the posterior is started. Our choice of 400 live 
points for the relatively well-behaved 7-dimensional distribution at hand 
is very conservati ve. In fact, as shown by the test problems examined in 
lFerozetalJ|2009l) , 1000 or 2000 active points are sufficient for the applica- 
tion of MultiNest even to, respectively, highly dimensional problems (with 
up to 30 parameters) or pathologically multimodal distributions (e.g. the 
egg-box likelihood presented in their Section 6.1). Moreover, if one is more 
interested in determining the mai'ginalized posterior distribution of the pa- 
rameters rather than in accurately calculating the value of the total evidence, 
a much smaller number of live points (of order 50) is shown to be already 
enough to obtain a reliable estimate. 



6 M. Barnabe et al. 




reconstructed source standard error 




X [arcsec] x [arcsec] x [arcsec] 

Figure 4. Best model lens image reconstruction for the system SDSS J0728. Top row: reconstructed source model; lo" uncertainty on the source pixels; 
significance of the reconstructed source. Bottom row: HST/ACS data showing the lens image after subtraction of the lens galaxy; lens image reconstruction; 
residuals. In the panels. North is up and East is to the left. 



of the lensed image map, where the residuals of the galaxy subtrac- 
tion are still appreciable. We use the best-fitting B-spline model of 
the lens galaxy as our data set for its surface brightness distribution, 
to avoid "contamination" from the background galaxy lensed im- 
ages, which is particularly bright and could bias the reconstructed 
model. A similar app roach had been follo wed in the analysis of lens 
system SDSSJ2321 jCzoske et all2008l) . 

The recovered non-linear parameters for the best recon- 
structed model, i.e. the maximum a posteriori model, are presented 
in Table[T] The uncertainties on the individual parameters are quan- 
tified by marginalizing, over each of them, the joint posterior dis- 
tribution (see Fig. The parameter values corresponding to the 
maximum of the one-dimensional marginalized posterior and the 
limits of the 99 per cent confidence interval are also listed in Ta- 
ble[Il 

We find for SDSS J0728 a logarithmic slope y' = l.OlT'^^l 
(errors indicate the 99% confidence level), very close to the so- 
called isothermal (i.e. p ~ l/r-) profile which appears to be a 
characterizing feature of early-type galaxies, and in general agree- 
ment with previous combined lensing and dynamics studies of 
the SLACS sample jKoopmans et al.l l2009l B09). For this spe- 
cific system — by using the 3-arcsec aperture averaged SDSS 
velocity dispersion me asure as the only kinematic constraint — 
iKoopmans et al.l | |2009|) determine a slightly different slope j' = 
1.85 ±0.10 (68% CL), which, however, we find here to be too shal- 
low to correctly reproduce the kinematic maps. One reason for this 
discrepancy is that the velocity dispersion derived from the SDSS is 
lower (typically of ~ 20 km s"') than the value obt ained from Keck 
spectr oscopy, as discussed in Sect. 12. 21 Moreover. iKoopmans et al.l 
( |2009|) adopt a simpler dynamical model, based on solving spher- 
ical Jeans equations, which might have difficulties in describing a 
system characterized by a significant flattening like SDSS J0728. 

The recovered lens strength is ao = 0323*^^^. The axial ra- 



tio of the total density distribution is found to be g = 0.602^° jj™, 
slightly flatter than the intrinsic axial ratio q^, = 0.688 of the lumi- 
nous distribution, calculated by deprojecting the observed isophotal 
axial ratio q^iD by making use of the best model value ; = 68.6^Q"g 
obtained for the inclination. Moreover, the position angle i?pa = 
67° 5 is found to be extremely close to the value inferred from the 
light distribution, indicating alignment between the dark and lumi- 
nous mass components in the inner regions of the galaxy. 

The marginalized posterior probability distributions of these 
parameters (with the exception of the position angle, which, as 
previously described, is kept fixed after the preliminary run), as 
well as those of the regularization hyper-parameters, obtained as 
described in Section [33] are shown in Fig.|6] We note that these 
represent the statistical errors on the considered power-law model, 
and do not take into account the potential systematic uncertainties 
due to issues in the generation of the dat a sets (e.g. the proced ure 
for galaxy subtraction, as pointed out by [Marshall et alj|2007l) or 
to incorrect model assumptions (as in the case of, e.g., non ax- 
ially symmetric density distribution or flattening that varies with 
radius). In the latter case, a quantitative estimate of the upper limits 
of the systematic errors can be obtained by looking at the findings 
of the lBamabe et all ( l2009bh ' crash-test' where the cauldron code 
is applied to a non-symmetric simulated galaxy. As mentioned in 
Sect. 13. ll the error on the logarithmic slope y' is less than 10 per 
cent even in this quite extreme case. Since real early-type galaxies 
are unlikely to depart from axisymmetry as drastically as this sim- 
ulated system, we expect systematic uncertainties to remain within 
a few per cent level. 

There is no evidence that any significant external shear effect 
needs to be included in order to model this system: when the exter- 
nal shear angle and strength are allowed to vary, the MAP model 
value of the latter is found to be negligibly small. This is consis- 
tent with the absence of massive structures near the lens galaxies 
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Figure 5. Best dynamical model for the galaxy SDSS J0728. Top row: observed surface brightness distribution, projected Hne-of-sight velocity and line-of- 
sight velocity dispersion. Middle row: coiTesponding reconstructed quantities for the best model. Bottom row: residuals. In the panels, North is up and East is 
to the left. 



found in iTreu et alj 120091) . We also explored the possibility of in- 
troducing a core radius in the density profile, without finding any 
improvement to justify the inclusion of this additional parameter. 

In Figure |4] and |5] we present — respectively for lensing and 
dynamics — the SDSSJ0728 data set, the reconstructed observ- 
ables obtained for the MAP model, and the corresponding resid- 
uals. This system displays an unusually structured lensed image, 
which we find to be produced by the presence of multiple compo- 
nents in the source plane. The reconstructed background source is 
shown in the top-left panel of Fig. |4j the top-middle and top-right 
panels show, respectively, the standard errors and the significance 
of the reconstructed source. 

The surface brightness map is reconstructed very accurately. 
Since the adopted data set is a B-spline model and therefore noise- 
less, the small residuals (at the 1 per cent level) are dominated by 
the discreteness effects of the TIC superposition, which determines 
the concentric ripples, while the structure aligned with the galaxy 



minor axis is caused by the toroidal shape of the building blocks. 
These undesired effects have been kept under control by increas- 
ing both the number of TIC employed in the dynamical model 
(of almost a factor of four, from Ne = 10 x Nl. = 10 elements 
adopted in the previous studies conducted with cauldron to 20 x 9 
elements) and the number of particles populating each TIC. An ad- 
ditional and more compelling motivation for this improvement of 
the TIC library, which further justifies the increased computational 
burden, is that it proves to be important for the reconstruction of the 
kinematic observables (particularly the velocity dispersion map) at 
larger distances from the galaxy major axis, i.e. that part of our data 
set which best allows us to probe the system under study beyond the 
elfective radius. 
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Figure 6. Marginalized posterior probability distributions of the power-law model parameters (inclination, lens strength, logarithmic slope, and axial ratio) 
and hyper-parameters, obtained from the nested sampling exploration of the posterior. The uncertainties on the parameters quoted in the text ai'e calculated by 
considering the interval around the peak which contains 99 per cent of the probability. 



3.5 Global dynamical quantities 

The recovered weighted stellar DF for the best combined model 
is presented in Fig. |7] as a map of the relative TIC weights over 
the integral space (£, L-) (the grid in the radial coordinate Rc is re- 
lated to a grid in energy as explained in BK07). This representation 
encodes, in a very compact way, much of the information on the 
dynamical structure of the galaxy that can be obtained under the 
assumptions of the adopted two-integral model. However, it is of- 
ten useful to distill such information into quantities that allow for a 
more straightforward physical interpretation. 

The global properties of the stellar velocity dispersion tensor 
are customarily encapsulated in the three anisotropy parameters 

p^l-'^, y^l-^ and 5^1- , (4) 

IT/;/; + n^^ 

where we indicate with 

Htt = J" pal d^x (5) 

the total unordered kinetic energy in the coordinate direction k 
and iTt(jr) is the velocity dispersion along the direction k (see 
ICappellari et alj2007l : lB"innev & Tremainell2008h . 

For SDSS J0728, we compute the integral of Eq. ^ within a 
cylinder of radius and height equal to R^, i.e. inside a region which 
is very well constrained by the data, finding a mild anisotropy 
5 = 0.08 ± 0.02, which falls within the typical range of values 
for early-type galaxies bot h in the local Universe ICappellari et al.l 
l2007l : lThonias et al.ll200^ and up to redshift z ~ 0.35 (see B09). 
Since we make use of a two-integral DF dynamical model, the ve- 
locity dispersion tensor is isotropic in the meridional plane (i.e. 
cr|(jc) = o":(x) for each x) and therefore p = by construction, 
and y is univocally linked to 6 by the relation y = 26/{S - 1). For 
the analyzed system, we have y = -0.18 ± 0.04. 

The importance of rotation with respect to random motions is 
among the most defining aspects of the dynamical structure of a 
stellar system. In order to explore how this property varies with the 



position in the meridional plane, we calculate the local ratio {v^)/d- 
between the mean rotation velocity around the z-axis and the mean 
velocity dispersion o"^ = {cr\ + a^ + al)l3 and we plot it in Fig.[8]up 
to one effective radius. The inner regions — within approximately 
1" — are dominated by random motions, while rotation becomes 
more important at large radii, a trend somewhat reminiscent of what 
is found, in B09, for the fast rotator SDSS J0959 (although it should 
be remembered that the {v^)l& map of the latter only extends up to 
Rc/2). Not surprisingly, therefore, the presence of large-scale or- 
dered motions is reflected also in the quite high value j- = 0.28^q 55] 
obtained for the intrinsic rotation parameter, which is a dimension- 
less proxy for the galaxy angular momentum (refer to B09 for the 
definition). 

3.6 Dark and luminous mass distribution 

The spherically averaged profile of the galaxy total mass corre- 
sponding to the best reconstructed model (solid black curve in 
Fig. |3.6| l can be calculated straightforwardly from Eq. ([TJ, while 
the radial profile of the luminous component is obtained from the 
recovered stellar DF. However, since within our method stars are 
treated as tracers of the total potential, the normalization of the lu- 
minous profile is not fixed, and must be constrained by means of an 
independent determination of the stellar mass-to-light ratio or by 
introducing additional assumptions. 

One particularly informative assumption, known as the "max- 
imum bulge" approach, consists in maximally reseating the lu- 
minous mass distribution without exceeding the total mass pro- 
file at any radius. This provides a consistent and robust way to 
assess a lower limit for the dark matter fraction in the analyzed 
system (under the hypothesis that the stellar mass-to-light ratio 
does not change too drastically with the position in the galaxy). 
By applying this procedure to SDSS J0728, we determine a value 
log(M,./MQ) = 11.50 for the total stellar mass. With this normal- 
ization of the luminous mass profile, one finds a dark matter frac- 
tion of 16 per cent of the total mass within the (three-dimensional) 
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Figure 7. Reconstruction of the weighted two-integral DF of the system 
SDSS J0728 obtained from the MAP model. The value of each pixel in the 
two-integral space represents the relative contribution of the corresponding 
TIC to the stellai' component of the modelled system. 



spherical radius r = R^jl, which rises to 28 per cent at r = 
— corresponding to about 6 kpc — and up to almost 40 per cent 
ai r ~ 10 kpc, which is approximately the outer limit of the re- 
gion over which we have direct information from the stellar kine- 
matic maps. This result is consistent with the findings of purely 
dynamical studies of t he inner regi ons of early-type galaxies in 
the local Universe (e .g. lOerhard et al . 2001; Cappellari et al. 2006; 
iThomas et al .l2007bh . It is also in good agreement with the conclu- 
sions of our previous analysis of six SLACS lens systems (where 
the same maximum bulge prescription was adopted), for which, 
however, the kinematic data set does not extend beyond R^, with 
the exception of galaxy SDSSJ0959 (see B09). The correspond- 
ing upper limit for the stellar mass-to-light ratio in the B band 
of SDSSJ0728 is M^,ILb = 3.14, which is in the lower end of 
the typical values r eported for slow-rotating elliptical galaxies in 
the local Uruverse jKronawitter et alj[200q ; lOerhard et"ai]|200ll ; 
iTruiillo et al]|2004l) . This is not too surprising, since the M-^ILg 
is expecte d to increase by a significant amount (a factor of 1.4 ac- 
cording to lTreu et al.ll2002l) between z = 0.2 and z = 0, simply due 
to passive evolution of stellar populations. Moreover, the evidence 
of rotation at large radii (cf . § 13.5b might indicate the presence of 
a disk component, typically characterized by a lower mass-to-light 
ratio. 

It is interesting to compare this maximum bulge upper limit 
with the stellar masses determined from stellar population analysis. 
By applying a novel Ba yesian stellar popul ation analysis code to 
multi-band imaging data, lAuger et"al] ( |2009|) determine for galaxy 
SDSSJ0728 — without including any priors from len sing — a 
value log(M*/Mo) = 1 1.44 + 0.12 for a lChabrieil ll2()03h IM F and 
a value log(Mj,/Mo) = 11.69 ± 0.12 for a lSalpeteil l ll955h IMF 
(quoted errors are 1-sigma). We note that, since the thermally puls- 
ing asymptotic giant branc h stars do not dominate the luminosity of 
ol d stellar popu l ations , the lBruzual & Chariot! ( l2003h models used 
by lAuger et al. I ( l2009l) should not be biased by ignoring them. In 
fact, the stellar masses for SLACS g alaxies are found to be consis- 
tent with the ma sses determined from lMarastorj | |2005|) models (see 
lTreuetal.ll2010l) . 

By using these stellar mass values to rescale the (spherically 
averaged) luminous profile, one obtains the red curves shown in 
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Figure 8. Map of the local {v^)/cr ratio between the mean rotation velocity 
around the j-axis and the mean velocity dispersion, plotted up to Re in the 
positive quadrant of the meridional plane 



Fig. 13.61 It can be clearly seen that the Salpeter IMF produces a 
stellar mass distribution which (unphysically) overshoots the total 
mass profile up to and beyond the effective radius. The Chabrier 
IMF, on the other hand, yields a physically acceptable luminous 
profile, which, interestingly, is also very close to the one deter- 
min ed through t he ma ximum bulge assumption. 

iTreu et ai] ( 1201 Oh argue - based on the mass determinations 
from lensing, dynamics and stellar populations synthesis mod- 
els for 56 SLACS systems - that early-type galaxies cannot have 
both universal IMF and dark matt er profiles. In fact, if a univer- 
sal N avarro, Frenk & White (NFW iNavarro, Frenk, & White! 19961 
1 19971) halo is assumed, the IMF shows a trend with velocity dis- 
persion: a Chabrier-like normalization is more appropriate for less 
massive systems with cr of the order of 200 km s"', while more 
massive galaxies are better described with a Salpeter-like IMF. The 
results of our detailed analysis of SDSS J0728 are consistent with 
this general picture. With a SDSS velocity dispersion cr = 214 ± 1 1 
km s"', this system definitely belongs to the lower mass end of the 
SLACS sample, and its halo profile is consistent with NFW (as dis- 
cussed in the next Section). Clearly, since this conclusion is based 
on only one system, a full analysis of the sample still needs to be 
done to further strengthen this tentative trend. 



3.7 Is the isothermal profile consistent with a NFW halo? 

The combined analysis of SDSS J0728 has provided us with the to- 
tal mass density profile of this galaxy within its inner regions. As 
discussed in the previous Section, and visualized in Fig. 13. 61 the lu- 
minous component alone cannot account for this mass profile over 
the entire radial range, even if its contribution is maximized (unless 
the stellar mass-to-light ratio changes with radius in a very fine- 
tuned manner). One, therefore, needs to invoke an additional mass 
component characterized by the specific profile that complements 
the luminous distribution. Interestingly, we find that the mass dis- 
tribution of this dark component can be consistent with a NFW halo 
profile. 

To show this, we attempt to describe the (spherically- 
averaged) total mass distribution of the galaxy as the sum of the 
luminous component, calculated as before from the recovered stel- 
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Figure 9. Spherically averaged mass distribution for lens galaxy SDSSJ0728. The solid black line shows the total mass profile obtained from the best 
reconstructed model, with the statistical uncertainty represented as a grey band. The sol id red line shows the luminous mass profile obtained from the recovered 
stellar DF and rescaled using the stellai' masses obtained from the I Auger et al.l )2009l) stellai' population analysis with Chabrier and Salpeter IMFs (left and 
light panel, respectively); the upper and lower eiTor bars quoted in that paper set the limits for the red-shaded regions. The vertical lines provide an indication 
of the region probed by the data set, showing the three-dimensional radius r which equals the Einstein radius (dashed line), the effective radius (dotted line) 
and the outermost boundary of the kinematic maps (dash-dotted line). 



lar DF, and a dark component modelled as a standard NFW halo, 
i.e. 



PnfwC'') ■ 



(6) 



where Pcnt is the critical density, the characteristic overdensity of 
the halo 6^ is a dimensionless parameter connected to the halo con- 
centration c by the relation 



6c 



200 



(7) 



3 [ln(l + c) - c/(l + c)] 

and i\ denotes the br eak radius. We adopt rj, = 25 kpc, based on the 
iGavazzi et alj ( |2007|) weak lensing analysis of the SLACS sample 
and the consideration that SDSS J0728 is slightly less massive than 
the average SLACS system (for which (ctsdss) - 250 km s"'). 

We then determine the (non-negative) normalization coeffi- 
cients for the dark and luminous distribution which allow to best 
reproduce (in a least-square sense) the total mass profile of the 
galaxy. Remarkably — despite the fact that the only two free pa- 
rameters here are the rescaling factors — the superposition of these 
two simple components proves to be enough to reproduce the total 
profile with great accuracy over the whole radial range covered by 
the observations, as shown in Fig.[TO] Moreover, the combined pro- 
file remains consistent with the one predicted by a nearly isother- 
mal density distribution even beyond 10 kpc and up to the break 
radius. 

The luminous profile obtained in this way almost coincides 
with the one determined by means of the maximum bulge approach, 
with log(My,/MQ) = 11.50, while the normalization for the dark 
halo profile translates into a concentration parameter c 11, from 
which one infers a virial radius ^200 = crt, ^ 280 kpc and a halo 
mass M200 - 3.1 x IO'-Mq. This is a mildly high value for the 



concentration when compared with the range c ~ 3 - 10 obtained 
from numerical simulations of rela xed dark matter haloes of corre- 
spond ing mass (see in particular lMaccio. Dutton. & van den BoschI 
I2OO8I) . However, concentrations higher than the theoretical predic- 
tions are found i n dynamic al studies of slow-rotating early-type 
galaxies (see e.g. iRomanow skv 2010, and references therein). We 
note that, in our case, lower concentrations are obtained by setting 
a larger value for the break radius, e.g. c ~ 10 for = 30 kpc, 
while the fit to the total mass profile becomes only slightly worse, 
with discrepancy of a few percent. 



4 CONCLUSIONS 

We have carried out a detailed study of the mass profile and dy- 
namical structure of the inner regions of the early-type lens galaxy 
SDSSJ0728, located at a redshift z = 0.21, using a composite 
data set consisting of HST/ ACS high-resolution imaging and two- 
dimensional kinematic maps constructed from long-slit spectro- 
scopic observations obtained with the Keck instrument LRIS (the 
slit — aligned with the galaxy major axis — has been positioned 
at three different heights along the minor axis, allowing to mimic 
integral-field spectroscopy). We have modelled the system by mak- 
ing use of the cauldron code for combined gravitational lensing 
and stellar dynamics analysis, which operates under the assump- 
tions of axial symmetry and two-integral stellar D F. With respect to 
sampl e studies of the SLACS lens galaxies such as lKoopmans et al] 
( l2009l) . the approach used here presents a number of improvements: 
it employs a self-consistent framework where the same total poten- 
tial is used for both lensing and dynamics; it allows one to construct 
genuine axisymmetric dynamical models (albeit restricted to two 
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Figure 10. Best-fit decomposition of the total mass profile of SDSS J0728 
(black curve) in the sum of a NFW dark matter halo with break radius r\y = 
25 kpc (blue curve) and the luminous profile obtained from the recovered 
stellar DF (green curve). The dotted red line which very closely matches 
the black curve over the whole radial interval is the profile obtained by 
adding up the luminous and dark components. The vertical dash-dotted line 
indicates the outermost boundary of the kinematic maps. 



integrals of motions); it extracts much more information from the 
data set, making use — in addition to the lensed image — of the 
lens galaxy surface brightness and velocity moments maps rather 
than being limited to a single measure of SDSS velocity disper- 
sion. This permits a much greater level of detail to be recovered 
when modelling the system. 

SDSS J0728 is the first galaxy for which it has been possible 
to conduct this kind of in-depth combined analysis by taking ad- 
vantage of a kinematic data set that extends well past the effective 
radius: the outermost pixels of the velocity moments maps probe 
a region up to a distance of about 3 arcsec from the center, corre- 
sponding to ~ 1.7 i?e- For comparison, in the sample of six SLACS 
lens galaxies examined in B09, the outermost boundary i?ki„ of the 
kinematic maps (obtained with VLT VIMOS integral-field spec- 
troscopy) is in the range 0.30 - 0.85 Rs, with the single exception 
of system SDSS J0959 for which R^i^ exceeds the effective radius 
of about 15 per cent. 

The main results of the analysis are summarized and discussed 
below: 

(i) The total mass density profile of the galaxy inner regions, up 
to a radius of ~ 1.7 R^, is found to be satisfactorily described by 
a simple axisymmetric distribution p oc 1 jni' , with a logarithmic 
slope y' = 2.08^q 552 (99 % CL). This is very close to the isothermal 
profile (i.e. p oc which seems to be an ubiquitous charac- 

teristic of early-type galaxies at least up to redshift z ~ 1, as re- 
vealed by a number of dyn amics, X-ray and combined lensing and 
dynamics studies (see e.g . Kronawitter et al. 2000I , Gerhard et al] 



I2OOIL iHumphrev & Buot3l2010l iKoopmans et al.ll2009l and refer 
ences therein). Moreover, weak lensing studies show that the total 
density profile remains consistent with being isothermal also well 



beyond the gal axy inner regions, up to radial distances of about 100 
effective radii ( IGavazzi et al.l2007h . 

The physical mechanisms which cause the total density distri- 
bution to be shaped into this particular structure, with little but 
non negligible intrinsic sc attering (<, 10 per cent in 7', see B09 
an d lKoopmans et al.ll2009h , are not well understood. The apparent 
"conspiracy" between the luminous and dark components to gen- 
erate a nearly isothermal combined profile appears even more sur- 
prising when considering that numerical simulations including gas 
dynamics find dark matter density profiles which depend both on 
the amo unt of baryons and on the details of the assembly processes 
(see e.g. lTissera et al.ll2009l : lAbadi et al.ll2009t) . 

Interestingly, despite the predicted complications, for 
SDSSJ0728 we find that it is possible to reproduce very ac- 
curately the (spherically averaged) total mass distribution by 
combining two very intuitive and simple building blocks: (1) the 
luminous mass profile obtained from the stellar DF, almost 
maximally rescaled and (2) a NFW dark matter halo (with a break 
radius r\, = 25 kpc and a concentration parameter c ~ 11). 

(ii) We find the total density distribution to be quite flattened 
within the probed region, with an axial ratio q = 0.60*^^1, which is 
flatter than the axial ratio q-^ of the luminous distribution (obtained 
by using the best model recovered inclination to deproject the two- 
dimensional isophotal axial ratio). This characteristic differentiates 
SDSS J0728 from the six lens galaxies studied in B09, for which 
q/qi, > 1. 

(iii) The system is characterized by a very mild anisotropy 
d = 0.08 ± 0.02. On examining the dynamical structure of the 
galaxy by means of the {v^}/d- map, one notices that the contri- 
bution of ordered motions becomes more important outside the in- 
ner regions, which determines the moderately high value /- = 0.28 
for the dimensionless angular momentum parameter. This result is 
obtained by integrating within a cylindrical region of radius and 
height equal to R^. If the integration is limited to R^/l, in order to 
allow a direct comparison with results of B09, one finds a lower 
j- = 0.18, fully consistent with the typical values obtained for 
the galaxies in that sample (with the exception of the clearly fast- 
rotating SDSSJ0959). 

(iv) Under the assumptions of maximum bulge and position- 
independent stellar mass-to-light ratio, we determine for the dark 
matter fraction a lower limit of 28 per cent within the spherical ra- 
dius r = R^. Within r ~ 10 kpc, i.e. the approximate extension 
of the area directly probed by the kinematic data, the contribution 
of the dark matter to the total cumulative mass is about 40 per 
cent, almost matching the luminous component. This is in agree- 
ment with the findings of dynamical studies of nearby ellipticals 
dOerhard et aLll200ll : Bappellari et alj|2006l : iThomas et ai]|2007al) 
as well as with the combined lensing and dynamics anal ysis of six 
SLACS galaxies at z = 0.08 - 0.33 jCzoske et al.ll2008l B09). In- 
terestingly, numerical simulations of early-type galaxy formation 
from cosmological initial conditions also predict a dark matter frac - 
tion of about 20 - 40 per cent within this radius jNaab et al .120071) . 

In alternative to the previous approach, we have also rescaled 
the luminous profile b y using the stellar masses calculated from 
the lAuger eTal] ( |2009|) stellar population analysis. The obtained 
luminous mass distribution is too high in the case of a Salpeter 
IMF, exceeding the total mass in places, while for a Chabrier IMF 
it remains lower than M,o,(r) and, moreover, close to the profile 
predicted under the maximum bulge hypothesis. This s uggests, 
in agreement with the conclusions of iTreu et ak I ( l2010l) . that a 
Chabrier functional form might be more suited to describe the IMF 
for less massive early-type galaxies such as SDSS J0728. 
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In order to further test if this description is correct, we plan to 
extend this study to a wider sample of SLACS lens galaxies, cov- 
ering a broad range of velocity dispersions cr w 200 - 350 km s"', 
for which two-dimensional kinematic data sets are available. 

In conclusion, despite being located at a redshift greater than 
0.2, the system SDSSJ0728 shows structural characteristics — 
namely nearly isothermal total density profile, dark matter frac- 
tion, anisotropy parameter 6, local ratio of ordered to random mo- 
tions — broadly consistent with what is observed in the nearby 
Universe for early-type galaxies of comparable luminosity and ve- 
locity dispersion (e.g . iGerhard et al.l l200ll iThomas et al.l 120073 . 
ICappellari et alj2007h . The upper limit for the B-band stellar mass- 
to-light ratio, obtained from the maximum bulge assumption, is 
also in line with the v alues determined for local ellipticals (e.g. 
iKronawitter et alfcOOOl) . once the ageing of the stellar populations 
is taken into account. This study, therefore, provides an indication 
that the density profile as well as the global dynamical structure of 
the inner regions of massive ellipticals did not undergo any dra- 
matic change or significant evolution across the last two billion 
years. 
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